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Abstract 

Flavour changing effects originating from the exchange of scalar particles in the pro- 
cesses e + e~ —>■ H x bs, H x bs, with H x = h° ,H° ,A°, are investigated in the Minimal Su- 
persymmetric Standard Model with non-minimal flavour violation at the one-loop level. 
The dominating SUSY-QCD contributions with squark-gluino loops are calculated and 
discussed. We consider the SUSY scenario with non-minimal flavour mixing in the down- 
type squark-mass matrix. The flavour- changing cross sections are derived, and we discuss 
the dependence on the MSSM parameters and the strength of flavour mixing. The values 
for the cross section can reach 10 -4 pb for the production of the heavy Higgs boson H° or 
A , and only 10~ 7 pb for the light Higgs boson h°. Non-decoupling behaviour occurs for 
both h°, H° production in the case of a common heavy SUSY mass scale. 
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1 Introduction 



Searching for Flavour Changing Neutral Currents (FCNC) is among the important tasks 
of the coming generation of high energy colliders [1] . FCNC processes are forbidden in the 
Standard Model (SM) at lowest order. These effects appear at the loop level and are of 
basic importance for testing the quantum structure of the SM. In the SM, however, the 
one-loop effects are small, suppressed by the GIM mechanism [2]. In models beyond the SM 
new non-standard particles appear in the loops, with significant contributions to flavour 
changing transitions [3]. Therefore, FCNC processes play an important role in searching 
new physics beyond the SM. 

Among various new physics models, supersymmetry (SUSY), especially the Minimal 
Super symmetric Standard Model (MSSM), is a favoured candidate [4]. The MSSM intro- 
duces two Higgs doublets to break the electroweak symmetry. After symmetry breaking, 
there are five physical Higgs bosons: two CP-even Higgs bosons (h°,H°), one CP-odd 
boson (A ) and two charged Higgs bosons {H^ [5]. The couplings of these Higgs particles 
may differ significantly from those of the SM. In fact, an important feature of SUSY models 
is that the fermion-Higgs couplings are no longer strictly proportional to the corresponding 
mass, as they are in the SM. In particular, the b-quark coupling to the neutral Higgs boson 
becomes enhanced for large tan/3 = v 2 /vi, the ratio of the two vacuum expectation val- 
ues [5] 0. Since there are five Higgs bosons in the MSSM, additional features not present 
in the SM may be useful to pin down differences between those models, manifested in 
different coupling strengths, decay widths and production cross sections. In the minimal 
flavour violation scenario (MFV), the only source of flavour violation is mediated by the 
CKM-Matrix [3]. In the general MSSM with non- minimal flavour violation (NMFV), new 
flavour changing (FC) effects emerge due to a possible misalignment between the squark 
and the quark sector. The FC interactions resulting from such a misalignment do not 
show up at the tree-level, but they can be generated at the one-loop level and could lead 
to relevant contributions to observables for specific regions of the MSSM parameters. The 
MSSM with NMFV is considered in this paper. 

The FCNC vertices in the SM and beyond have been extensively examined in the 
literature, and the results promise FCNC to provide a fertile ground for testing the SM and 
probing new physics. These analysis include rare decays of 5-meson systems [6,7], Z-boson 
decays [8-10], and top decays [11]. It is known that the SM predictions for the top quark 
FCNC processes are far below the detectable level and that the MSSM can enhance them 
by several orders, making them potentially accessible at future collider experiments [12]. 
In particular, the branching ratios for Z —>■ bs decays are of the order of 10 -8 in the SM [8], 
of the same order in SUSY with t-c mixing and can reach 10~ 6 in SUSY with b-s mixing [9], 
both last rates being dominated by the SUSY radiative effects from squark-gluino loops. 
FC effects that can be induced by squark-gluino loops in Higgs-boson decays have been 
investigated as well, and BR(H X — > bs) ~ 10~ 4 — 10~ 3 and BR(H X — > tc) ~ 10~ 4 have been 
found for selected regions of the MSSM parameter space and the flavour mixing parameters 

Explicitly they are proportional to 1/ cos /?, and to [— sin a , cos a, sin /?] for [h° , H° , A°) respectively. 
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in the MSSM squark-mass matrices [13,14]. The electroweak corrections are sub dominant, 
at least one order of magnitude smaller than the SUSY-QCD ones, but can give interference 
effects [15]. The large rates found for the SUSY contributions to the Higgs partial decay 
widths into bs and sb, as well as to the effective FC Higgs couplings to quarks, are indeed 
quite encouraging [13-19]. Furthermore, effects of non-minimal flavour violation on the 
MSSM Higgs boson masses and the electroweak precision observables, m w and sin#e/7, 
has been investigated in [20]. 

FC effects in the SM in electron-positron collisions have been under study in several 
papers. In particular, the cross sections of e + e~ — > tc and e + e _ — > bs processes in the SM 
are analyzed in [21]. The cross section of the process with a bs final state was found to 
be larger than the one with a tc final state, but only of the order of 10 -3 fb, which is too 
small to be observed at present and future colliders. 

In this paper we investigate FC effects via associated bottom-strange and Higgs boson 
production in the MSSM. More concretely, we study the effects of squark-gluino loops in 
the higher-order processes 

e+e" -> H x bs + H x bs (H x = h° , H° , A ) , (1.1) 

and discuss the size of the induced cross sections at the one-loop level, with virtual second 
and third generation squarks, and their dependence on the MSSM parameters and the 
flavour mixing strength. 

The paper is organized as follows: In section [2] we present a brief outline on squark 
mixing in the MSSM with NMFV. Section [3] gives an overview over the various classes of 
squark-gluino diagrams at the one-loop level and their impact on the total cross section, 
including both resonant diagrams (subsection 13. 2p and non-resonant contributions (subsec- 
tion [33]). We discuss the dependence on the MSSM parameters in section |U Additionally, 
the non-decoupling behaviour of the SUSY contributions in the large sparticle-mass limit 
is presented in section |5l and a summary is given in section [6j 

2 Non-minimal squark mixing 

Flavour changing phenomena in SUSY models can emanate from a mixing between different 
generations of squarks through the soft breaking terms in the Lagrangian of the squark 
sector. This non-minimal flavour mixing scenario, where squark mixing results from a 
misalignment between the quark and the squark mass matrices, is the most general case 
of the MSSM. We assume that the flavour changing squark mixing is significant only in 
transitions between third and second generation squarks. It is known that there are strong 
experimental bounds on the squark mixing involving the first generation, resulting from 
data on K° - K° and D° - D° mixing [3]. 

In general, the flavour violating quantities arise from non-diagonal entries in the bilinear 
soft breaking matrices M~ , Mj~ and Mg, that appear in the mass matrices in the up- 
squark and down-squark sectors as well as from non-diagonal (flavour changing) entries in 
the trilinear soft breaking matrices A u and Ad- In our analysis, FC effects are generated 
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by the one-loop exchange of 6-s-admixture states and gluinos, g. When the squark mass 
matrix is diagonalized, FC gluino- quark- squark couplings can be derived for the mass 
eigenstates and FC effects via squark-gluino loops are induced. It is known that some of 
the flavour mixing parameters in the b — s sector are severely constrained by the b — > sj 
decays [3,17,22,23], but in particular the ones referring to the LL- and RR-mixing of the 
SUSY partners of the left-handed and right-handed quarks, respectively, are not definitely 
excluded [3]. 

We assume the simplest case of mixing, where the only non-zero off-diagonal squark 
squared mass entries in the squark-sector stand for s^l and c^l mixing. Expressing the 
squark mass matrix in terms of the chirality eigenstates {sl, sr, &l, &r} and {cl, cr, £l, £r}, 
respectively, and considering only LL-Mixing, i.e. mixing between the left-chiral compo- 
nents of the squarks leads to the following mass matrix mediating between the chiral 
eigenstates and the down-squark admixtures [13,15], 
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(2.3) 



where M^, Mr a are the corresponding bilinear soft SUSY breaking entries and m q the 
quark matrices with q G {b, s,c, t}. T| is the third component of the weak isospin, Q q 
denotes the charge of the corresponding quark. The trilinear soft SUSY breaking matrices 
A q are assumed to be diagonal, with entries A Sjb . \i is the mass parameter of the Higgs 
boson sector and tan (3 is the ratio of the vacuum expectation values in this sector, mz is 
the Z boson mass, and sw = sm9w contains the electroweak mixing angle 9w- A = Xll 
denotes the parameter characterizing the flavour mixing strength in the mixing of the LL- 
components of the second and third generation squarks 0. Thus, we have only one free 
parameter determining the strength of FC. 



2 Thc flavour mixing parameter A can be identified with (Sll)23 in the usual notation of the mass- 
insertion approximation [3]. 
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3 The loop-induced cross sections 



In this section we discuss the various contributions leading to the production of a Higgs 
boson associated with a flavor-nondiagonal quark pair, bs or sb. As far as the final states 
with bs and sb are experimentally not distinguished, the result for a(e + e~ — > H x b s+H x s b) 
is obtained by summing the two individual cross sections. 

We have done the computation using FeynArts, FeynCalc and FormCalc [24]. Feynman 
rules of MSSM vertices with FC effects are recently implemented in FeynArts (extending 
the previous MSSM model file). The Higgs boson masses, the Higgs boson decays and the 
masses of the MSSM-spectrum have been computed using the FeynHiggs2.1beta Fortran 
code [25], including all FC effects. As the MSSM Higgs boson couplings to down- type 
fermions receive large quantum corrections enhanced by tan/5, these corrections have to 
be resummed to all orders in perturbation theory with the help of the effective Lagrangian 
formalism [26]. These leading threshold corrections to the bottom mass are taken into 
account in our computation in terms of the resummed version. 



3.1 Kinematics 

The cross section for the production process 

e'ipa) + e + { Ph ) -> 6(pi) + s{p 2 ) + H x {p 3 ) 
is given by the expression (neglecting the electron mass) 



a = 



2s(2tt) e 
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i=i 



d 3 Pi 
2rf 



5\p a + p b -Y,P*) l-^l 2 ' 



(3.1) 



(3.2) 



involving an integration over the 5-dimensional phase space and \M.\ 2 , the spin-averaged 
matrix-element squared of the process. The phase space integral can be expressed in terms 
of two energy and two angular variables, chosen in the overall CMS with s = (p a + p b ) 2 , 
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where 



cr = \fs~ — p\, t = a 2 — |pl| 2 , m± = m 2 ± m 3 



9 is the angle between the momenta p a and pi, and r] the angle between the two planes 
spanned by pi,p 2 and pi,p a - 
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(d) 

Figure 1: (a), (b), (c) Generic diagrams for resonant contributions to a(e + e~ —>■ <pb s+<p s b) 
(0 = h° , H° , A ), (d) Gluino-squark loop contributions to Z — > b s and <fi — > bs. 



3.2 Resonant contributions 

Due to the structure of the couplings of Higgs bosons to another Higgs boson and/or to a 
Z-boson in the MSSM, the qqH production processes receive large resonating contributions 
not present in the SM. The associated Higgs boson production with a heavy quark pair, 
tt and bb, in e + e~ collisions at high energy have been investigated and large resonant 
contributions that imply rates of a few fb have been found [27] for the 66-pair. Here we 
consider the resonating contributions to the associated Higgs boson production with a 
quark pair 6s in e + e~ collisions. 

The structure of the corresponding Feynman diagrams is shown in Fig. [TJ other di- 
agrams are suppressed by the small electron mass. In the kinematically allowed regions 
there are two types of resonating intermediate states for H x = h°, H° production, involv- 
ing Z and A resonances, e + e~ — > ZH X with Z —* bs, and e + e _ — > A°H X with A — > bs 
(Fig. []Ja,b)). There are other resonances channels for A production, involving h° and H° 
resonances; namely e + e~ — > A°H X with H x — > bs (Fig. d](c)). The result for the H x = h° 
case (h° resonance) is several orders of magnitude smaller than the H° resonance case 
(see below), therefore we do not include this case in the discussion. These processes are 
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sensitive to the MSSM couplings of the Z boson to the Higgs boson H x and to the A°H X 
pair, normalized to the standard ZZHsm coupling given by 

9zA°h° = cos(/5 - a) , g ZZh o = sin(/? - a) , 

9 z ao ho = sin(/3 - a) , g ZZH o = cos((3 - a) . (3.5) 

For the computation, we treat the intermediate particle in the narrow-width approxi- 
mation, i.e. multiplying the two-particle production cross section for on-shell particles with 
the corresponding branching ratios, 

a(e + e- -f H x bs + H x sb) ~ cx(e + e~ -> H x Z) ■ BR(Z -> bs + sb) , 
a(e + e~ ^ H x bs + H x sb) ~ a(e + e~ — > A ) • 6 s + s 6) , 

a(e + e- ^ H x bs + H x sb) ~ a{e + e - ^ H x A ) ■ BR{H X ^ bs + sb) , (3.6) 

which yields a good approximation owing to the small total decay widths of the interme- 
diate resonances. The flavor-changing Z, A and .£P decays are loop-induced, with the 
gluino-squark strong contributions shown in Fig. QJd) as the dominating source, yielding 
the branching ratios of Y — Z, A , H x , 

B *(y^, + .i ) = r < y -"> + r(^'*> , (3.7) 

1 tot ^ J 

where r tot (V) is the total width in each case. 

Although numerical results and a discussion of the parameter dependence will be given 
in the next section, we want to illustrate here the size of the various contributions, choosing 
rriA = 250 GeV, tan/3 = 35, and the SUSY parameters of (14.11) . For the value A = 0.6 of 
the mixing parameter, one obtains the following branching ratios, 

BR(A° -> bs + sb) ~ 10~ 2 , BR(H° -^bs + sb) ~ 10~ 2 , 

BR(h° ^bs + sb) ~ 10" 4 , BR(Z^ bs + sb) ~ 10~ n . (3.8) 

Here the threshold corrections on the 6-quark mass are considered. This implies that the 
results for Z — > bs are different from the ones given in [9]. For mj = 5 GeV we can 
reproduce the previous results. 

The couplings given in (13.51) can be suppressed if either sin(/3 — a) or cos(/5 — a) is 
very small. In particular for large tan/3 and not too low tua,, as it is the case for the 
choice of input parameters (14.11) . the ZA°h° and ZZH° couplings are suppressed and the 
couplings for ZA°H° and ZZh° are large. Hence, the cross sections for e + e~ — > H° A 
and e + e~ — > h° Z are about 2-3 orders of magnitude larger than for e + e~ — > h° A and 
e + e~ — > H° Z, at y/s = 500 GeV amounting to 

a(e + e- -> H° A ) ~ 10~ 2 pb , a(e + e" -> h° Z) ~ 10~ 2 - 10" 1 pb, 
a(e + e~ -> ~ 10~ 5 - 10~ 4 pb, a(e + e - -> #° Z) ~ 10~ 5 - 10~ 3 pb . (3.9) 
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Figure 3: Resonant contributions to a(e + e — > H x bs + H x sb) as a function of A at 
yfe = 500 GeV, 1000 GeV, 1500 GeV, for h° (left), H° (right). 

Therefore, we can conclude that the Z-resonance contribution for the H x bs final state 
(H x = h° , H° , A ) are very small, at most of the order of 10~ 12 pb, and the /i°-resonance 
contributions are of the order of 10 -8 pb, whereas the contributions with heavy Higgs 
bosons (A , H°) intermediate states can reach 10 -4 pb for the cross section. 

In Fig. [2] we show the behaviour of the cross section with the center of mass energy, 
based on the dominant contributions with the A intermediate states, Fig. [T](b), for both 
h° (left panel) and H° (right panel). Since the results for the contributions with H° 
intermediate states, Fig. [He), are very similar, we do not include explicit numerical results 
for this case. The behaviour of the cross section with both the center of mass energy and 
the flavour parameter A for the case of A intermediate states are straightforwardly applied 
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to the H° case. If not stated differently, the value of A is fixed to be A = 0.6. The variation 
of the cross section with A is contained in Fig. [3j shown for three different values of the 
energy, y/s = 500 , 1000 , 1500 GeV. The effects increase with A for both kind of Higgs 
bosons and vanish at A = 0, as to be expected. In the H° case, due to kinematical reasons, 
the cross section at y/s = 500 GeV is smaller than at 1 TeV, leading to an inversion of the 
order of the curves. In particular for the light Higgs boson h°, the cross section remains 
rather small, although it is effectively for a two-particle process with subsequent decay. 
The resonating behaviour is largely canceled by either a small Zh°A° coupling or a tiny 
branching ratio of the Z boson. Therefore, also non-resonant contributions have to be 
considered for a reliable estimate of the total production rate. 

3.3 Non-resonant contributions 




Figure 4: Generic diagrams for the squark-gluino one-loop contributions to e + e — > H x b s+ 
H x sb (H x = h°,H°,A°). 
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Figure 5: One-loop contributions to cr(e + e~ — > H x bs + H x sb) as a function of \fs (GeV). 
Total cross sections, main and box contributions, and corrections with Higgs radiation off 
the s-lines (h° = left panel, H° = right panel). Here A = 0.6. 

Besides the two-particle-like processes discussed in the previous subsection, there are 
non-resonating contributions for H x bs production, which are genuine 2 — > 3 processes. The 
set of diagrams under consideration is shown in Fig. HI There are three different types of 
relevant topologies: self-energy, triangle, and box diagrams. For the self-energy diagrams 
one has to distinguish between two different cases of Higgs radiation, one with radiation 
from Mines and the other one with radiation from s-lines. Because of the small Yukawa 
couplings for the s case, only Higgs radiation from 6-quarks will be of significance. 

Here we outline the size of the contributions of each class of diagrams to the cross 
section for illustrational purposes, based on the same set of model parameters as in sub- 
section I3.2L The non-resonant e + e~ — > Z — ► bs contributions with Higgs radiation from 
b, s-lines (Fig. IHj,k)) are suppressed with respect to the other contributions, being of the 
order of 10~ 12 — 10~ 14 pb, therefore they are not included in the analysis. 

In Fig. [5] the one-loop contributions from the diagrams in Fig. H] are depicted as a 
function of the center of mass energy y/s. Since the results for a(e + e~ — > H° bs + H° s b) 
and a(e + e~ — > A b s + A s b) are very similar, we do not include explicit numerical results 
for A production. The discussion for the H° case can straightforwardly be applied to the 
A case along this paper. Besides the total cross sections (solid line) we list in addition the 
contributions from three different topologies independently: main contributions (short- 
dashed line), defined as the ones resulting from triangle diagrams and from self-energy 
diagrams with Higgs radiation off the b- fines, box contributions (long-dashed line), and the 
s-rad contributions, defined as the self-energy contributions with Higgs radiation off the 
s-lines. These last contributions are presented here for illustrational purposes, they are of 
the order of roughly 10 -11 pb for both Higgs bosons. In the case of the heavy Higgs bosons, 
the box contributions are smaller than in the h° case, being of the order of 10~ 10 — 10~ n pb, 
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therefore they appear together with the contributions with Higgs radiation on the s-lines 
near to zero. For the case of h° production, there are several peaks corresponding to 
squark-mass thresholds in the loop integrals. These peaks are less distinctive in the H° 
case because they appear only in the box contributions, which are smaller for H°. 

One can see from Fig. [5] that for both Higgs bosons the non- resonant contributions are 
roughly two orders of magnitude smaller than the resonant contributions analyzed before. 
The dependence of the cross sections on the mixing parameter A is illustrated in Fig. O at 
1 TeV. As before, the cross sections increase with A, being exactly zero for A = 0. 

4 Dependence on the MSSM parameters. 

The SUSY parameter set needed to determine the input for e + e~ — > H x bs + H x sb 
(H x = h°, H°) consists of the quantities tua, tan (3, fi, Mg, Mo, A, where Mo is a common 
value for the soft-breaking squark mass parameters, M = Mg g = Mfj, c ^ = Mg, / Sj6 \, and 
A denotes the trilinear parameters which are chosen to be equal, A t = Ab = A c = A s = A; 
M g is the mass of the gluino. These six parameters will be varied over a broad range, 
subject to the requirements that all the squark masses be heavier than 150 GeV [28] 
and Mh > 114.4 GeV [29]. Similarly, in view of the present experimental bounds on 
the chargino mass [28], we consider only values above 90 GeV. The flavour mixing 
parameter, A, is constrained by the lower squark-mass bounds. For those parameters that 
are kept fixed in the various figures the following default set has been chosen, 

/i = 1000 GeV , M = 500 GeV , A = 800 GeV , 
m A = 250 GeV , M g = 800 GeV , tan = 35 , (4.1) 
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which is in accordance with experimental bounds for the decay b — > 37 [30] , as checked with 
the help of the code micrOMEGAs [31], based on minimal flavour violation calculations 
at leading order [32] and some contributions beyond leading order that are important for 
high values of tan/3 [33], as well as in accordance with the result for the muonic (g — 2). 
The FCNC effects in the rare decay b — > 57 may impose severe constraints on MSSM 
parameters that reduce the maximally allowed value for BR(h° — > bs) [14]. Furthermore, 
starting from (14. ip . A > 0.8 implies unallowed values for the squark masses, i.e. Mg < 150 
GeV. If not stated differently, A = 0.6 is a default value. 

In the following we explore the dependence of the cross section on the MSSM parameters 
for both resonant and non-resonant diagrams, illustrated in Figs. [7] to Figs. [TQl therein, 
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Figure 8: Cross section as a function as a function of \i (GeV) at \J~s = 500 , 1000 , 1500 
GeV for: (a) e + e~ —> A°H X — > H x bs + H x s b, (b) non-resonant e + e~ — * H x bs + -fP s 6. 



the upper panels correspond to the resonant contributions and the lower panels show the 
non-resonant contributions. 

First, we discuss the dependence on tan/3, the parameter that appears inside the cou- 
plings as well as in the squark-mass matrices and thus determines the squark-mass splitting 
and the b-mass corrections. Fig. [7| contains the cross section cr(e + e~ — ► H x bs + H x sb) 
versus tan (3 for three values of the center mass energy, yfs = 500 , 1000 , 1500 GeV. In 
all cases, the cross sections increasing with tan j3. For the h° boson (left panel) the in- 
crease with tan f3 is slower due to a superposition of two different effects in the involved 
subprocesses: the reinforcement in the decay subprocess from the coupling enhancement 
is partly compensated by decreasing in cos {(3 — a) in the production cross section of the 
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Figure 9: Cross section as a function of Mg (GeV) at y/s = 500 , 1000 , 1500 GeV for: (a) 
e + e~ — > A°H X — > H x bs + H x s b, (b) non-resonant e + e~ — > H x bs + H x s b. 



2^2 process, which causes the almost linear shape of the tan (3 dependence. Values of 
tan/3 < 10 are not included in these plots because they imply Higgs boson mass values 
lower than 114.4 GeV. 

The dependence on the \x parameter for the h° (left) and H° case (right) is shown in 
Fig. [HJ The shaded regions in these figures correspond to the regions excluded by the LEP 
bounds on the chargino mass, < 90 GeV. One can see that in all cases the results 
for the cross sections are not symmetric under // — > This is essentially due to the 
inclusion of the threshold corrections to the bottom quark mass, where the sign of /x plays 
an important role. We found that the corrections grow with the /i parameter. 

Next we study the behaviour with respect to the other MSSM parameters. It turns out 
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Figure 10: Cross section as a function of M (GeV) at y/s = 500 , 1000 , 1500 GeV for: (a) 
e + e~ — > A°H X — > H x bs + H x s b, (b) non-resonant e + e _ — > H x bs + H x s b. 



that the cross sections are nearly independent of the common trilinear parameters A. In 
Fig. [9] we show the dependence on the gluino mass, M- g . One can see from this figure that 
for all cases the FCNC effects vanish for heavy gluinos, i.e. the gluino decouples when only 
the gluino mass is varied. It is also interesting to see that the maximum is reached not for 
the smallest values but for intermediate values of the gluino mass. 

The dependence of the cross sections on the common soft SUSY breaking squark mass 
parameter M GeV is displayed in Fig. [TD1 Regions below M = 400 GeV are not dropped 
because they would lead to excluded values for the down-squark masses. The FC cross 



15 



sections decrease with M , showing the independent decoupling behaviour of M . 

In summary, the study of the dependence on the MSSM parameters, has confirmed 
the maximum size of the cross sections of the order of 1CT 7 pb for h° and O(10~ 4 ) pb for 
the H° case and demonstrates the individual decoupling behaviour of the virtual SUSY 
particles. 

5 Non-decoupling behaviour of heavy SUSY particles 

Here we analyze the non-decoupling behaviour of squarks and gluinos in SUSY-QCD con- 
tributions to our FCNC e + e" processes that result in a non-vanishing cross section also 
for very heavy SUSY particles. Non-decoupling behaviour of SUSY particles also occurs 
in the SUSY-QCD and SUSY-EW contributions to FCNC Higgs boson decays [13,15]. 
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Figure 11: Non-decoupling behaviour of the cross section from resonant contributions 
(upper panels), and of a(e + e~ — > H x —+ bs + sb) from non- resonant contributions (lower 
panels). Ms, in GeV, is defined in (15.11) . 



16 



In general, a non-decoupling behaviour originates from a compensation of the mass sup- 
pression induced through the heavy-particle propagators, by mass terms in the interaction 
vertices. Non-decoupling behaviour has previously been analyzed in detail in the case of 
flavour-preserving MSSM Higgs boson decays [34-37]. The non-decoupling contributions 
to effective FC Higgs Yukawa couplings have also been studied in the effective-Lagrangian 
approach for the quark sector [18,19] and the lepton sector [38]. 

A condition for non- decoupling of SUSY particles is a common scale for all SUSY mass 
parameters. For the simplest assumption, the values of the SUSY masses, where all the 
soft breaking squark mass parameters, collectively denoted by M , the \i parameter, the 
trilinear parameters, collectively denoted by A, and the gluino mass Mg, are chosen to be 
of the same size and much larger than the electroweak scale Mew, 

M s = M = M g = fj, = A > M EW ■ (5.1) 

The numerical results of the analysis are shown in Fig. [11] for both h°, H° production 
(for A = 0.6). One can see that the cross sections tend to a non- vanishing value for large 
values of Ms, a behaviour characteristic of non-decoupling. While Ms is displayed only up 
to 2 or 3 TeV, we have checked that the cross sections remain stable even for very large 
values of Ms- 

6 Conclusions 

We have considered the production of neutral MSSM Higgs bosons in association with bs 
and sb quark pairs. We analyzed the flavour changing effects emanating from squark-gluino 
one-loop contributions, which are the dominant ones in a scenario with non-minimal flavor 
mixing in the squark sector. Both resonant contributions with an intermediate on-shell A 
boson and non-resonant contributions have been investigated. The highest cross sections 
for e+e" -> h°bs + h°sb are of 0(1(T 7 ) pb and of O(10" 4 ) pb for e+e" -> H°bs + H°sb 
(and similar for A ) which will be too low to be seen at a future linear collider. 

One feature of the SUSY-QCD contributions of general basic interest is their non- 
decoupling behaviour for large values of the SUSY particle masses. The flavour changing 
cross sections, when all soft breaking mass parameters, \x and the trilinear parameters, are 
of the same order of magnitude and much larger than the electroweak scale, do not vanish 
but tend to a constant value. 
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